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EFFECTS OF UNSTEADY FREE STREAM VELOCITY AND FREE STREAM 
TURBULENCE ON STAGNATION POINT HEAT TRANSFER 


SUMMARY 

An analysis is presented to study the combined effects of 
transient free stream velocity and free stream turbulence on heat 
transfer at a stagnation point over a cylinder situated in a 
crossflow. An eddy diffusivlty model has been formulated and the 
governing momentum and energy equations are integrated by means of 
the steepest descent method. The numerical results for the wall 
shear stress and heat transfer rate are correlated by a turbulence 
parameter. It has been found that the wall friction and heat transfer 
rate increase with increasing free stream turbulence intensity. 
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NOMENCLATURE 


Cf = skin friction 

D = diameter of cylinder 

f = frequency of disturbance 

^ = dimensionless frequency of oscillations (fD/2Lto) 

k = thermal conductivity 

K = constant (3.631 IWD) 

Nu = Nusselt number 

Pr = Prandtl number 

Pr + = turbulent Prandtl number 

q = heat flux rate 

Re = Reynolds number (UoD/v) 

T = temperature 

Tu = turbulence Intensity 

u,v = velocity components 

U m = reference velocity 

x = coordinate In streamwlse direction 

y = coordinate normal to the surface 

p = density 

a = thermal diffusivity 

T) = dimensionless distance 

t = dimensionless time 

0 = dimensionless temperature 

a w = wall shear stress 

ft = dimensionless frequency (irfD^U^) 

Subscripts 

s = steady state 

w = wal I conditions 


00 


ambient conditions 
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INTRODUCTION 

Turbomachinery blades provide some of the most challenging pr6b!ems 
of fluid mechanics. The presence of secondary effects in the boundary 
layers over rotating blades results in the flow chracteri sties that 
are different from the classical two-dimensional boundary layer theory. 

The extent of three dimensionality depends on the angular velocity, 
the flow coefficient, space chord ratio, aspect ratio, stagger angle 
of the blade, etc. Turbine rotor blades rotate at speeds up to 15,000 rpm 
In the wakes of several stationary vanes. The precise role of the free 
stream velocity fluctuations superimposed on the complex flows that occur 
In turbomachinery blades is not yet completely understood. 

There exists a need for a systematic and detailed theoretical and 
experimental study of the effects of the upstream flow disturbances on 
the flow and heat transfer characteristics of a gas turbine blade. The 
flow disturbances at the exit of the stationary vanes are pulsatile in 
character and so the velocity and temperature distribution as well as the 
local skin friction and heat transfer coefficients become transient. In 
multistage turbomachinery, the wake of a blade following the first row 
of blades always evolves under the influence of the turbulence present 
In the wake of the previous row of blades. 

Many modern turbine blade sections, even at the high Reynolds 
numbers at which they operate, manifest considerable areas over which 
the boundary layer remains laminar. The superimposition of mainstream 
velocity f lucturations associated with artificially induced mainstream 
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turbulence is known to affect dramatically the laminar fluid flow and 
heat transfer characteristics. The life of gas turbine blading depends 
to a large extent on the local heat transfer rates to the blade surface 
from the free stream fluid. A study of the combined effects of the 
mainstream turbulence and the unsteady mainstream velocity on the flow 
and heat transfer from a gas turbine blade is therefore of considerable 
practical importance. 

Sutera et al. Cl, 23 simulated the free stream turbulence by a distri- 
buted vorticity in the free stream. Their model assumed a vorticity 
component in the oncoming flow to be unidirectional and oriented so that 
the vortex lines were susceptible for stretching. They found that ampli- 
fication by stretching of vorticity of sufficiently large scale can occur. 
As a result, such vorticity, with small intensity in a flow can appear 
near the boundary layer with an amplified intensity. They also found 
that the thermal boundary layer is apparently much more sensitive to 
the induced effects than is the velocity boundary layer. Further eluci- 
dation of the mechanism by which the presence of the vortices affects 
the heat transfer has been provided by Kestin and Wood [[3,43. Smith 
and Kuethe C53 postulated an expression for the eddy diffusivity induced 
by the free stream turbulence. Galloway [[63 proposed a model on the 
assumption that the heat transfer in the presence of free stream tur- 
bulence is enhanced by Goertler vortices. Traci and Wilcox [[73 pre- 
sented a two-equation model of turbulence. Early theoretical investi- 
gators [8,93 predicted the effect of the free stream oscillations on 



5 


the mean boundary layer heat transfer or shear stress. These studies 
resulted in a good qualitative understanding of the flow mechanisms 
Involved, but they do not provide the capability for quantitative 
predictions. 

J unkhan and Serovy CIO] reported experimental data concerning the 
effects of free stream turbulence intensity on the flow and heat 
transfer from an isothermal flat plate with a favorable pressure 
gradient. Their data indicates that there is no effect of free 
stream turbulence intensity on heat transfer through a laminar 
boundary layer with zero pressure gradient. 

There exist several practical applications in which wakes evolve 
in the presence of free stream turbulence. One such situation exists 
in single or multistage turbomachinery, in which the wake of a blade 
following the first row of blades always evolves under the influence 
of the turbulence present in the wake of the previous roiHr of blades. 

An extensive search of literature on general two- and three-dimensional 
wakes is reported by Raj Of], Pal and Raj Cl2] have investigated the 
effect of free stream turbulence on the characteristics of a turbulent 
wake developed from the trailing edge of a thin and smooth flat plate 
both analytically and experimentally. They found that the free stream 
turbulence increases the wake recovery and growth rates. 

The present study was undertaken in order to investigate the combined 
effects of free stream turbulence and time dependent free stream 
velocity on the unsteady flow and heat transfer characteristics of a 
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stagnation point on a circular cylinder. This basic study should 
enable understanding of one important aspect of the complex flow and 
heat transfer phenomena in the real turbomachinery environment to a 
first degree of approximation. 

Gorla, Jankowski and Textor E I 3H investigated the time-mean 
characteri st ics of the laminar boundary layer near an axisymmetric 
stagnation point when the velocity of the oncoming flow relative to 
the body oscillates. They presented different solutions for the small 
and high values of the reduced frequency parameter. Gorla [14] 
proposed a model for the eddy diffusivity induced by free stream 
turbulence intensity and integral length scale. The eddy diffusivity 
model was applied to the stagnation point of a ciruclar cylinder situated 
in a uniform crossflow in the presence of free stream turbulence. A 
numerical solution of the governing momentum and energy equations 
yielded results for the skin friction coefficient and the Nusselt 
number. The effects of free stream turbulence intensity and integral 
length scale on the local steady state heat transfer rate from a 
circular cylinder at various angles around its periphery measured from 
the stagnation point were studied by Gorla and Nemeth [15]. Correlations 
were provided by Gorla [16] for the experimental data reported by Base 
et al. [17] concerning the heat transfer at the stagnation point under 
the combined influence of turbulence intensity and free stream velocity 
oscillations. A detailed analysis was made of the combined effects of 
transient free stream velocity and turbulence intensity on the unsteady 
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stagnation point flow characteristics by Gorla Cl8] while these effects 
on unsteady stagnation point heat transfer are reported recently by 
Gorla D9U. 

Computer time for this project was provided by the Computer Center 
of the Cleveland State University. A detailed literature review was 
performed by Mr. Robert Platteter. 
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GOVERNING EQUATIONS 


The flow development and the coordinate system are shown in 
Figure 1. Assuming an incompressible flow with constant properties 
and negligible dissipation, the governing equations, within the 
framework of the boundary layer approximation may be written as 
3u 3v n 

Mass: -k + 37 = 0 

i 3u . 3u . 3u 3U. . .. 3Un , 1 3o 
Momentum: at + *J^+''^ = _ 3f+U e - § f+-J7 


Energy: + u 21 + v — = - — 

3t u 3x v 3y pC p 3y 


(1) 

( 2 ) 

(3) 


In the above equations, x and y represent the distances along the 
streamwise direction and normal direction; u and v the velocity components 
in x and y directions; U e (x,t) the external velocity; t the time; 
p the density of the fluid; T the temperature, Cp the specific heat of 
the fluid; o the shear stress and q the heat flux. 

The boundary conditions for the wall and for the outer edge of the 
boundary layer are: 

y=0: u=v=0 and T = T w 
y -*■ »; u = U e (x,t) and T-*Teo 


( 4 ) 
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To complete the formulation of the problem, initial conditions must 
be specified in the (t,y) plane for x = 0 and in the (x,y) plane for 
t = 0. Consideration will be given here for a flow in which at time 
t = 0, the flow field is given by steady-state conditions. At time 
0» +he external velocity U e (x,t) begins to deviate from the steady 
state velocity Uq(x) 

U e (x,t> = U 0 (x) *G(t) (5) 

In the above equation, G(t) represents unsteadiness of the free stream. 

It is assumed that G(t) = 1 for t < 0. 

The shear stress and the heat flux in equations (2) and (3) may be 
written as 


where 


° - P <v + % ) 57 

q = - pC p (« +p ^>. §1 

£ m = momentum eddy diffusivity 


Pr-f-= turbulent Prandtl number 


(6) 

(7) 


It has been assumed that the conventionally used value of the turbulent 

Prandtl number Pr + = 0.9 is applicable to the present problem. 

Proceeding with the analysis, we define a stream function ip such that 

, dip 

u - 7 ^ and v = - ma Y be verified that the continuity equation is 

automatically satisfied. We further define 
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( 8 ) 


Substitution of expressions in (8) into the momentum and energy equations 
yields the following nondimensiona I equations for the velocity and 
temperature fields: 

• • • 

(1ts)f ,,, +(s»+f-| 2 )f ,f -{ ^(f’-D+g' +(f r ) 2 -1> (9) 


1 m + Pr ' s } - ( n6 

Pr ci+ T F7 ) 9 J - 2G^ 


where 

S = ( e m /v) . 

The transformed boundary conditions are 

f ' ( t, 0) = f ( t,0) - 0, f'Cx, 00 ) = 1 
6(t,0) = 1 , 6(x,“) = 0 


( 11 ) 


In the above equations, primes indicate differentiation with respect 


to q and the dot with respect to x. 
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EDDY DIFFUSIVITY FORMULATION 

Upon substituting ^ = 0 and G=1 into equations (9) and (10), 
one obtains governing equations for the velocity and temperature fields 
under steady state conditions. The momentum and energy equations for 
the steady-estate problem are given by 

[(i+s)f"]'+ff«'+ci-(f) 2 : = o (i2) 

+ p^ ),0 '3 ,+ f 0f = ° <i3) 

The appropriate boundary conditions are given by 
f(0) = f'(0) = 0 and f'(~) = 1 

6(0) = 1 and 0(») = 0 (14) 

The steady-state friction factor and Nusselt number N u may be calculated 
after solving equations (12) and (13) and may be written as 

C f = (W>- '3-838 "5) 

Nu = - 1.906 Re^ 0 T (O) (16) 


where 


Re 


^ UooD ^ 
v v ; 


In equations (12-16), it must be noted that f = f(q) and 0 = 0(q) only. 

Equations (12) and (13) are solved on a computer using the fourth- 
order, Runge-Kutta numerical procedure. The double precision arithmetic 
was used in all the computations. The selection of the integration step 
size depends upon the desired level of accuracy. After some computational 
trials, a step size of q = 0.001 was chosen. 
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In the present paper, we restrict ourselves to the case of homogene- 
ous and isotropic free-stream turbulence. We have assumed the following 
expression for the eddy dlffusivity in the present work. 

~ = s = si (n+n 2 ) U7) 

where s^ = constant proportional to the turbulence intensity 
and n = dimensionless distance normal to the surface. 

A review of the literature suggests that there are several uncertain- 
ties in the published experimental data including tunnel blockage, vari- 
ation of physical properties, effect of turbulence scale, and inaccuracies 
in the measurement of turbulent intensity. Restricting the discussion to 
the stagnation point, it has been found convenient to use the single 

i 

correlation parameter (Tu*Re 2 ). This parameter has also been suggested 
by Smith and Kuethe C5H on the basis of a semi -emp i r i ca I theory. 

Figures 2 and 3 show the results for the steady-state friction 
factor and heat transfer rate, respectively, as functions of the 
correlation parameter (Tu*Rei). The constant s^ in the present eddy 
diffusivity model has been adjusted so as to have a best fit of the 
present numerical predictions with experimental data available in the 
! Iterature. It has been found that the following expression for s-| 
makes the predictions of the present analysis to be in good agreement 
with the published data over the entire range of (Tu-Re^): 

St = 0.018 (Tu-Re*) (18) 

The results from the present analysis are seen to agree with the 
experimental data within 15 percent. 
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SOLUTION 


The governing momentum and energy equations (9) and (10) for the 
unsteady problem are integrated by the steepest descent method. 

Considering the momentum equation first, we let 

f(T,n) = Z a n ( T )- 2? (19) 

n=2 

Substituting the expression for f(T,n) into equation (9) and collecting 
coefficients of like powers of n, we find that 

G 

33 = -( s-| 32 + 1 + q2 ) 

• • • 

a 4 = 1^ + f §2 ' a 2 + 2s 1 2a 2 + 2s 1 ^ 1 + §2 5 
a 5 = 3s 1 a 4 + a 2 + ^ 

a 6 = ~^ s ] a 5 + ^^2 a 4 + +2a2a3 (20) 


etc. 


The unknown coefficient a 2 ^x) is determined through the use of 
boundary and initial conditions. When once a 2 <x) is determined, the 
remaining features of the flow field such as the transient velocity 
profiles may be readily evaluated. The time-dependent wall shear stress 
and friction factor may be written as 


°W ^ ( 3y^x,o,t 


= y - kx GC-^) 2 .a 2 (x) 




( 21 ) 
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We now define an integrating factor: e- c ^ T » ri ^ 

n i n A 

where I = / (y^-) * ( 5 *+f -^1 dTi 

Momentum equation then becomes upon integration 
f » * = e 1 • a(x,Ti) 

n x * * 

where a = a 2 (x) + / (f’-D+^t (f') 2 -0 dp 

In order to arrive at the boundary condition at the outer edge of the 
boundary layer, momentum equation is integrated from q = 0 to 00 . 

Thus 

1 = / f »» dr) = / e -I(T ' n) ‘ dr 1 

0 0 


To evaluate the Integral In equation (23), the series for f(x,r|) 

to be inverted. We now define a new variable y such that 

,, _ t / » i- ^ r q 2 Sin 3 Si 2 n 5 -| r ln /2si — si 2 \ 2 

Y - I(T,q) + qZ Lj — + io“J- siq-( 5 } n 


oo 

= E Cn n 
n=0 


n+3 


where C 0 = ff “ f s i : 


Cl = 22. - s . 1 a ? + 

4! 8 4 

C 2=ft-^ + |s 1 3 -fs 1 4 + ^- 

= <s, 4 -2s, 3 H2s,-2| z )] 


(22) 


(23) 

I s 

(24) 
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Inversion of series for y in equation (24) is given by 


r| = E A, 




m=0 (m+1) 


* Y 


(m+1 )/3 


(25) 


Upon applying Cauchy’s residue theorem, we get 

fl 1 . - ( m+ 1 ) /3 . 

A m = 2?] / Y dn 

where i = J~^~\ . The integrations are carried out three times in the y plane 
to dispose of the fractional powers of y once in the q plane. We now have 


where 


1 = / f" dq = / e • i|>(n,T) 
o o 


* = 2 fr*n" 

n=0 n ' 


dn 

dy 


dy 


(26) 


B o = a 2 

. 

P / J+Gy / 1 , G N 

B 1 = -<gT-> (1+s 1 “ 4 q 2^ ~ s 1 a 2 

B 2 = (1+S! - 4§Z>'(§2 a 2 + | ) +2si(1 + §2 )(1+ Si - 4§7> 

+2a 2 ^4§Z ~ S 1 + S'] 2 ) 

B 3 = 2(1+ Sl -•j 7 )(||J+i|+a2 2 )-3si(1+s 1 - j?) (§2 *2 + ^> 

G , , . . G „ , G 


-6(1 +g2)d+ Sl - Jz^Jz-si +S! 2 ) 

+6a 2 (^- - s 1 3) 


(27) 


etc. 
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From equation (26) we recall that 

V =y -2 ^ 3 ’ E Y m ^ 3 

m=0 


where 


d o = ^ ' c o' 


d 1 -i 


1/3 

2/3 


1 -2/3 2 Ci 

c 0 * (B 1 3 B 0 c Q ) 


d 2 = J [-B 0 C 2 C 0 - 2 +B 0 C) 2 C 0 - 5 -BiCiC 0 -2 + f- C 0 ~ ] J 
d3 = 3 B 0 c 3 c 0 " 7/3 + f- b 0 c 1 c 2 c 0 " 10/3 

-^ C ^ BoCo - 13 / 3 -^,^- 7 / 3 

+ BjC^Co- 1073 -! B 2 C 1 C 0 “ 7/3 + ^ C 0 ' 4/3 : 


etc. 


Equation (26) will now be Integrated In terms of the complete 
gamma functions to yield 


it 

3n 


CO 

(T,“) =1 

m=0 


dm 


r( 



(28) 


Upon applying Euler's summation procedure to the above series, we have 
1 =jf • r( j)’ d o(T> +^T(f)-d 1 (T) 

+ -jj • d 2 (T) + -j^-r(j) • d 3 (T) ( 29 ) 

After substituting expressions for dg, d^, 62 and 6 ^ into the above equation, 
we obtain a first-order ordinary differential equation for 32 (t). Equation 
(29) has been integrated numerically by means of a fourth-order, Runge-Kutta 
procedure. A discussion of the results will be presented in the next 
section. 
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It should be noted that the transient velocity profiles may be 
determined when once a 2 (x) is determined. The velocity profiles are 
given by the equation 

It = E dj T (J±t <30) 

3n J=0 Y 3 


where r, is the incomplete gamma function, 
by equation (21 ) . 

In order to integrate the unsteady energy 

0° p) 

0(T,n) = s bn ^ ) ‘ n 

n=0 n ‘ 


The friction factor is given 


equation (10), we set 


(31 ) 


where 


b 2 = -Pr s^ bi 


b 3 = PrQ 1 


-2Prsibi + Pr 2 si 2 bi 


(32) 


etc. 

In equation (32), we have 

q (33) 

y 1 Pr + 

Proceeding as before, the energy equation reduces to the following form: 

9GPr + Co 4/5 + dbj_ _ 1 . 9GPr t Op / # {A Q r(i) 

Pr s^ T( 4/3) dT 3 Pr si T ( 4/3) 

t AiT(|) t A 2 > • b-| ( T> 


(34) 
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where 


C A = C— _ 19.1 - s .1 + filG-i 

0 L 6 3Pr+ 6G^ J 


'0 

A o 


6 3Pr + 
1/3 


ITT-) 

°0 


A 1 = - 


24Cl . 5/3 

— (J-) 

36 ^0 


a 2 = “ 


C 2 (C 1) 2 

2 + 6 2 

(C 0 ) 


( c 0 ) 


p, _ £1 Ql32 
0/1 " 
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Ho = SA. _ 9m 4 . . ,28^, G . (Qi 4 -2(?3 S1 

~ 120 30 + 5 Pr + 2G^ ) + 5Pr+ 


(35) 


The time-dependent wall heat flux and Nusselt number may be 
written as 

% ■ -Kf<£> 

x,o,t 


= -Kfdw-Tco) SW 0»(O,x) 

V 

N u = -D/KG/v 

when once bj(T) Is calculated, the temperature profile as well as the 
wall heat transfer rate may be deduced. Equation (34) was numerically 
Integrated by means of the fourth-order, Runge-Kutta procedure. 



DISCUSSION 


It must be noted that the analysis Is performed for predicting 
the transient response behavior at the stagnation point due to the 
combined effects of an arbitrary time dependent free stream velocity 
and homogeneous, isotropic free stream turbulence. No assumption has 
been made for G(t) that represents the unsteadiness of the free 
stream. In fact, the present analysis is valid for any arbitrarily 
varying continuous G(t). There exist in the I iterature some experi- 
mental data and some predictive methods to analyze separately the 
wall shear factor at a stagnation point under steady-state conditions 
In the presence of free-stream turbulence and the friction factor due 
to transient free stream velocity with no free stream turbulence. 
Neither experimental nor analytical work has been reported to investi- 
gate the combined effects of unsteady free stream velocity and free 
stream turbulence at a stagnation point. 

In order to illustrate the application of the present analysis 
to investigate the title problem and also to assess its accuracy, 
we have chosen two different cases of time-dependent, free stream 
velocity. 

The first case corresponds to G(t) = (1+t) with no free stream 
turbu I ence present. Physically, this implies that at time t=0 the 
inviscid flow begins to change with constant positive temporal 
acceleration. Yang H20U investigated this problem be approximate 
integral methods. The results for a 2 (T) from the present analysis 
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are compared with those of Yang £20] in Fig. 4. The quasi-steady 
solution given by Yang incorrectly predicts a discontinuity in the 
wall shear at t = 0. The results from the present analysis start at 
a 2 (0) = 1 .23259. After reaching a maximum value of a 2 initially as t 
increases, a 2 decreases with increasing t towards the quasi-steady 
solution. The present analysis is seen to yield physically more 
realistic results than Yang's integral analysis for x< 1, during 
which time the flow is experiencing its greatest unsteadiness. The 
agreement between the present results and those of Yang is within 
2 percent for t > 1.0. 

The second case corresponds to G(r) = 1+A sin(fix). This represents 
a stagnation point flow In which the free stream oscillates with 
amplitude A and dimensionless frequency 12 about a mean velocity. 

Lighthiil £8] studied this problem using a perturbation method and 
integral approach. L i ghth ill’s analysis was valid for fi ->■<». Ishigaki £9] 
studied the same problem retaining second-order terms in the perturbation 
analysis. The results from the present analysis for the shear phase 
advance are compared with those of Ishigaki and Lighthiil in Fig. 5. 

The agreement between the present results and those of Ishigaki was 
seen to be within 5 percent. Also presented in Fig. 5 are the present 
results for wall hat flux phase lag. L i ghth ill’s asymptotic result as 

is also in Fig. 5. 

Figures 6 and 7 illustrate variation of the shear stress function 
a 2 (t) under the combined effect of transient free stream velocity and 
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free stream turbulence. For the sake of illustration, G(t) = 1 + A sJn(ftr) 

only was chosen. It is seen that any existing transients die out within 

1 oscillation of the free stream. After about 1 oscillation, the shear 

stress function a 2 ^T) assumes a periodic shape. The amplitude of the 

1 

oscillations Is a strong function of the turbufence parameter (Tu-Re 2 ). 

It must be noted that ^( t ) increases substantially from the steady-state 
solution for wall shear stress. Physically, a 2 (x) corresponds to the wa! 
shear stress at given instant of time as given by equation (15). It may 
be observed that even after all starting transients have died out the 
wall shear does not fluctuate with equal amplitude about the steady-state 
value. 

Figure 8 illustrates the variation of the normalized wall heat 
transfer under the combined effect of transient free stream velocity 
and free stream turbulence. We have chosen G(t) = 1 + sin(ftT) for this 
purpose. It is seen that any existing transients die out within two 
oscillations of the free stream. After about two oscillations, the 
wall heat flux assumes a periodic shape. The amplitude of the oscilla- 
tions was observed to be a strong function of the turbulence parameter 
(Tu-Rei). It may be observed that even after all starting transients 
have died out the wall heat flux does not fluctuate with equal amplitude 
about the steady-state value. 

It must be noted that the dimensionless frequency parameter ft 
may be written as (iTfR/UU) where f is the frequency of disturbance, 

R the cylinder radius and U® the free stream velocity. The grouping 
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(fR/Uj is the Strouhal number. Ishigaki C2lj derived from his theo- 
retical analysis that the parameter J^Tu^ is a control ling parameter 
to describe heat transfer in turbulent oscillating flows. Bayley 
and Priddy C22] reported heat transfer measurements around a gas 
turbine blade by controlling the frequency and amplitude of turbulence. 
The Nusselt number data was correlated by then in terms of a parameter 
(fC/U„) 2 . Tu^*Re where C was the chord length of the blade. For the 
title problem including both unsteady free stream velocity and free 

stream turbulence, a 2 (T) could not correlate well with the parameter 

1 

(Tu’Re 2 ). It may be recalled that the steady-state friction factor 
data has been successfully correlated by the parameter (Tu*Re^) as 
shown in Fig. 2. Utilizing the analogy between heat and momentum 
transfer mechanisms, it was hypothesized that (^Tu 2 Re) may be a suit- 
able correlation parameter. 

Figure 9 illustrates a comparison between the present predictions 
and the experimental results of Base et al. [J7H for the dimensionless 
heat transfer rate versus the correlating parameter f 2 Tu 2 Re. Although 
the present results underestimate the Nusselt number at high values of 
f 2 Tu Re, the agreement between the two seems to be more encouraging at 
lower values of the correlating parameter. 
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CONCLUDING REMARKS 

The method of approach used in this paper to the analysis of the 
combined effects of free stream transient velocity and free stream 
turbulence at a stagnation point reduced the governing equations to 
first-order ordinary differential equations. The method is capable 
of yielding results for any arbitrary time-dependent nature of the free 
stream velocity. The solutions are not restricted by small perturbation 
assumptions. A successful formulation has been performed to obtain an 
expression for the eddy diffusivity induced by the free stream turbulence 
intensity. The solution of the governing unsteady momentum and energy 
equations with the proposed eddy diffusivity model yielded predictions 
for the skin friction coefficient and heat transfer rate. A correlation 
parameter is suggested to correlate the friction factor and heat transfer 
results for the title problem. It has been found that the wall friction 
and heat transfer rate increase with increasing free stream turbulence 
intensity. In the case of flows involving unsteady free stream velocity, 
the friction factor and heat transfer rate increase with increasing 
values of the reduced frequency of oscillations. 
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FIGURE 2. Effect of Turbulence Intensity on Friction Factor 
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FIGURE 3. Effect of Turbulence Intensity on Heat Transfer Rate 
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FIGURE 4. Sheer Stress Function Versus r for G(r) = (1+r). 
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FIGURE 5. Sheer Phase Advance and Wall Heat Flux 
Phase Lag for G(t)= 1 + A sin (Qr) 
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FIGURE 8 Normalized Heat Transfer Versus r 
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Experimental Data of Base et al 



FIGURE 9 

Nusselt Number Versus the Turbulence Parameter 
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